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Abstract: We work out and review prospects for future quarkonium physics 

Q_i' at HERA. We focus on the impact of color-octet contributions and discuss 

^ ■ how measurements at HERA can be used to test the picture of quarkonium 

■ production as developed in the context of the NRQCD factorization approach. 

<N ' 

^ 1 Introduction 

o 
in 

The production of heavy quarkonium states in high energy collisions provides an impor- 
tant tool to study the interplay between perturbative and nonperturbative QCD dynam- 
ics. Quarkonium production in deep inelastic scattering and photon-proton collisions at 
HERA has been analysed at some length in the context of the previous HERA workshops 
(see e.g. ref.jl]]). However, most of the previous studies have been carried out within 
the color-singlet model (CSM) || or the color-evaporation model ||. Only recently, a 



X 



rigorous theoretical framework for treating quarkonium production and decays has been 
developed in ref. Qj. The factorization approach is based on the use of non-relativistic QCD 
(NRQCD) H to separate the short- distance parts from the long-distance matrix elements 
and explicitly takes into account the complete structure of the quarkonium Fock space. 
This formalism implies that so-called color-octet processes, in which the heavy-quark an- 
tiquark pair is produced at short distances in a color-octet state and subsequently evolves 
nonperturbatively into a physical quarkonium, should contribute to the cross section. Ac- 
cording to the factorization formalism, the inclusive cross section for the production of 
a quarkonium state H can be expressed as a sum of terms, each of which factors into a 
short- distance coefficient and a long-distance matrix element: 

da(ep ^H + X)=Y, Mep -> QQ M + X) (O h [n]) . (1) 

n 

Here, da denotes the short-distance cross section for producing an on-shell QQ-pair in a 
color, spin and angular-momentum state labelled by n. The NRQCD vev matrix elements 
(O h [n]) give the probability for a QQ-pa.ii in the state n to form the quarkonium state 
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H . The relative importance of the various terms in (JJ) can be estimated by using NRQCD 
velocity scaling rules ||. For v — > (v being the average velocity of the heavy quark in 
the quarkonium rest frame) each of the NRQCD matrix elements scales with a definite 
power of v and the general expression can be organized into an expansion in powers 
of v 2 . 

It has recently been argued in refs.JF], [| that quarkonium production in hadronic 
collisions at the Tevatron can be accounted for by including color-octet processes and 
by adjusting the unknown long-distance color-octet matrix elements to fit the data. In 
order to establish the phenomenological significance of the color-octet mechanism it is 
however necessary to identify color-octet contributions in different production processes. [] 
In this report we will focus our discussion on the prospects of extracting information on 
the color-octet processes from the measurements of inelastic quarkonium production at 
HERA. Elastic/diffractive mechanisms will be discussed elsewhere [IIJ. We shall briefly 



review the impact of color-octet contributions and higher-order QCD corrections on the 
cross section for J/ip photoproduction and compare the theoretical predictions with recent 
experimental data. In this context we will also comment on the possibility of measuring 
the gluon distribution in the proton from J/ip photoproduction. The high-statistics data to 
be expected in the future at HERA will allow for a detailed comparison of the theoretical 
predictions with experimental data not only for J/ip photoproduction, but also for various 
other channels and final states, like photoproduction of ip', T and x states, associated 
J/ip + 7 production, fragmentation and resolved-photon contributions as well as deep 
inelastic J/ip production. We will discuss how these reactions can be used to constrain the 
color-octet matrix elements and test the picture of quarkonium production as developed 
in the context of the NRQCD factorization approach. 



2 J/0 photoproduction 

Quarkonium production in high energy ep collisions at HERA is dominated by photo- 
production events where the electron is scattered by a small angle producing photons of 
almost zero virtuality. The measurements at HERA provide information on the dynamics 
of quarkonium photoproduction in a wide kinematical region, 30 GeV <^ y/s 1P <, 200 GeV, 
corresponding to initial photon energies in a fixed-target experiment of 450 GeV <^ E 1 
20, 000 GeV. The production of J/ip particles in photon-proton collisions proceeds predom- 
inantly through photon-gluon fusion. Elastic/diffractive mechanisms can be eliminated by 
measuring the J/ip energy spectrum, described by the scaling variable z — p • kip I p • fc^, 
with p, being the momenta of the proton and J/ip, 7 particles, respectively. In the 

1 Quarkonium production via color-octet states has also been studied in e + e~ annihilation and Z 
decays, for hadronic collisions in the energy range of fixed-target experiments and in B decays ||. 



2 



proton rest frame, z is the ratio of the J/0 to 7 energy, z = E^/E^. For elastic/diffractive 
events z is close to one; a clean sample of inelastic events can be obtained in the range 
z^0.9. 

For J/0 production and at leading order in v 2 , the general factorization formula ([I]) 
reduces to the standard expression of the color-singlet model |§ . The short- distance cross 
section is given by the subprocess 

l + g^cc[ 3 S u l}+g (2) 

with cc in a color-singlet state (denoted by 1), zero relative velocity, and spin/angular- 
momentum quantum numbers 2S+l Lj = 3 Si. Relativistic corrections due to the motion of 
the charm quarks in the J/0 bound state enhance the large- z region, but can be neglected 
in the inelastic domain [|TTJ. The calculation of the higher-order perturbative QCD cor- 
rections to the short- distance process fl2|) has been performed in refs.flT^, [13|| . Inclusion 
of the NLO corrections reduces the scale dependence of the theoretical prediction and 
increases the color-singlet cross section by more than 50%, depending in detail on the 



photon-proton energy and the choice of parameters [13 



Color-octet configurations are produced at leading order in a s through the 2 



parton processes [14, 15, II 



1 + g - cc[ l S ,%] 

1 + g - cc[ 3 P ,2,8]. (3) 

Due to kinematical constraints, the leading color-octet terms will only contribute to the 
upper endpoint of the J/0 energy spectrum, z ~ 1 and pr ~ 0, px being the J/0 transverse 



momentum. It has, however, been argued in ref. [17| that sizable higher-twist effects are 
expected to contribute in the region p T <^ 1 GeV, which cause the breakdown of the 
factorization formula ([!]). Moreover, diffractive production mechanisms which cannot be 
calculated within perturbative QCD might contaminate the region z Ri 1 and make it 
difficult to extract precise information on the color-octet processes. 

It is therefore more appropriate to study J/0 photoproduction in the inelastic region 
z < 0.9 and px > 1 GeV where no diffractive channels contribute and where the gen- 
eral factorization formula (0) and perturbative QCD calculations should be applicable. 
Color-octet configurations which contribute to inelastic J/0 photoproduction are produced 
through the subprocesses [^D, [U| 



7 + # -> cc[ 1 S ,8]+g 
j + g — cc[ 3 S 1 ,8]+g 

j + g — cc[ 3 P ,i,2,8] +g. (4) 
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Light-quark initiated contributions are strongly suppressed at HERA energies and can 
safely be neglected. Adopting NRQCD matrix elements consistent with those extracted 
from the fits to prompt J/ip data at the Tevatron || (see Table P one finds that color- 
octet and color-singlet contributions to the inelastic cross section are predicted to be 





1.16 GeV 3 


m 3 f 3 


(OW[ l S ,8]) 


10- 2 GeV 3 


3 7 




lO- 2 GeV 3 


3 7 

m c v 




10-2 GeV 3 


1 7 



Table 1: Values of the NRQCD matrix elements used in the numerical 
analysis, with the velocity and mass scaling, v is the velocity of the 
heavy quark in the quarkonium rest frame. For charmonium it holds 
v 2 ~ 0.25. 

of comparable size [14], [L6|. However, taking into account the uncertainty due to the 
value of the charm quark mass and the strong coupling, the significance of color-octet 
contributions can not easily be deduced from the analysis of the absolute J/ip production 
rates. A distinctive signal for color-octet processes should, however, be visible in the 
J/i[) energy distribution da/dz shown in fig.[l| ||. Since the shape of the distribution is 
insensitive to higher-order QCD corrections or to the uncertainty induced by the choice 
for m c and a a , the analysis of the J/ip energy spectrum da/dz should provide a clean test 
for the underlying production mechanism. From fig.|l] we can conclude that the shape 
predicted by the color-octet contributions is not supported by the experimental data and 
that the J/ifj energy spectrum is adequately accounted for by the color-singlet channel. 
This is however not to be considered as a failure of the factorization approach, since also 
other analyses (see for instance the last item of ref. |J) have pointed out that the fits to 
the Tevatron data may have returned slightly too large values for the matrix elements. 
With higher statistics data it will be possible to extract more detailed information on the 
color-octet matrix elements, in particular from the analysis of the J/ip energy distribution 
in the inelastic region. 

The impact of higher-order QCD corrections on total cross sections and differential dis- 
tributions has been studied thoroughly for the color-singlet channel in ref. JT3|. A detailed 
analysis of the spectra in the high energy range at HERA shows that the perturbative 
calculation is not well-behaved in the limit pt — > 0, where pr is the transverse momentum 
of the J/0. No reliable prediction can be made in this singular boundary region without 
resummation of large logarithmic corrections caused by multiple gluon emission. If the 
small pt region is excluded from the analysis, experimental results on differential distri- 
butions and total cross sections are well accounted for by the color-singlet channel alone, 
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Figure 1: Color-singlet and color-octet contributions to the J/ip energy 
distribution da / dz at the photon-proton centre of mass energy y/s~ p = 
100 GeV integrated in the range > 1 GeV [14]. Experimental data 
from O, 111. 



including NLO QCD corrections, see e.g. fig.|2|. However, since the average momentum 
fraction of the partons is shifted to larger values when excluding the small-pr region, the 
sensitivity of the prediction to the small-x behaviour of the gluon distribution is not very 
distinctive. A detailed analysis reveals that the size of the QCD corrections increases 
when adopting parton densities with flatter gluons. The sensitivity to different gluon 
distributions is thus reduced in next-to-leading order as compared to the leading-order 
result, in particular when choosing a small charm mass and a large value for the strong 
coupling. Parametrizations with extremely flat gluons like MRS(DO') |21] are clearly dis- 
favoured by the recent HERA measurements of the proton structure function and do 
not allow for a reliable prediction in the high energy region. For the parameters adopted 
in fig.0, the MRS(DO') distribution leads to next-to-leading order results not very different 
from those obtained with the MRS (A') parametrization. The corresponding i^-factors are 
however uncomfortably large, K ~ 4, casting doubts on the reliability of the perturbative 
expansion as obtained by using flat gluon distributions. If parton distributions with steep 
gluon densities are adopted, the next-to-leading order cross section is well-behaved and 
gives an adequate description of the experimental data, as demonstrated in fig.0. 
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Figure 2: NLO color-singlet prediction for the total inelastic J/ip pho- 
toproduction cross section as a function of the photon-proton energy for 



different parametrizations 2C] of the parton distribution in the proton 



13 1 . Experimental data from [18, 



The cross section for the production of ijj' particles has been measured by several photo- 
and hadroproduction experiments [23|] to be suppressed by a factor ~ 0.25 compared to J/ip 
production. This result is consistent with naive estimates obtained within the color-singlet 
model where one assumes that the effective charm masses in the short distance amplitudes 
scale like the corresponding ip' and J/ip masses ||. Within the NRQCD factorization 
approach it is however conceptually preferred to express the short-distance coefficients in 
terms of the charm quark mass rather than the quarkonium mass, yielding a(ip') /a(J/ip) ~ 
0.5 for color-singlet dominated production channels. Relativistic corrections and color- 
octet contributions are expected to affect the ratio o \ip') / 'a (J /fy>), but no quantitative 
prediction can be made with the present experimental and theoretical information. 

Bottomonium production is a particular interesting subject to be studied at HERA. 
The larger value of the bottom quark mass makes the perturbative QCD predictions of 
the short-distance cross section more reliable than for charm production. Moreover, the 
derivation of the factorization formula eq. ([!]) relies on the fact that the momentum scales 
which govern the bound state dynamics are well separated: (rriQV 2 ) 2 <C {rnqv) 2 <C mq. 
This assumption is reasonably good for charmonium (where v 2 ~ 0.3) but very good 
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for bottomonium (where v 2 ~ 0.1). Thus, the theoretical predictions for bottomonium 
production should be much more reliable than those for charmonium. The production 
rates for T bound states are, however, suppressed, compared with J/ip states, by a factor 
of about 300 at HERA, a consequence of the smaller bottom electric charge and the phase 
space reduction by the large h mass. 



3 J/i/j production via fragmentation 

At sufficiently large transverse momentum p?, quarkonia production is dominated by 
fragmentation, the production of a parton with large px which subsequently decays into 
the quarkonium state and other partons While the fragmentation contributions are 
of higher order in a s compared to direct quarkonia production in fusion processes, they 
are enhanced by powers p T /m 2 c and can thus overtake the fusion contribution at pt ^ m c- 
It has in fact been argued that quarkonium production at large pt in hadronic collisions 
at the Tevatron can be accounted for by including gluon fragmentation into color-octet 
states 0| . 

The fragmentation contribution to the differential cross section for producing a quarko- 
nium state H at large pt can be written in the factorized form 

d<7 ( 7 + p^H + X)=J2 r dC d<x(7 + P -> i + X) D^ H {Q (5) 

where da is the differential cross section for producing a parton of type i with momentum 
Pt/C- The fragmentation function gives the probability that a jet initiated by 

parton i contains a hadron H carrying a fraction £ of the parton momentum. According 
to the NRQCD factorization formalism, fragmentation functions for charmonium have the 
general form 

A-^(C, a0 = £ d^ cs[n]x (c, (o H N) (6) 

n 

analogous to eq.(|l|). The function di^ c5 [ n ]x gives the probability for the parton i to form 
a jet that includes a cc pair in the state labelled by n. It can be calculated at an initial 



scale /i ~ m c as a perturbative expansion in a s {m c ) |24|, g5[ and be evolved up to higher 



scales fi ~ pt by using the Altarelli-Parisi evolution equations. 

A quantitative analysis of the fragmentation contributions to J/ip photoproduction 



at HERA has been performed in ref.||26|| (see also ref. ||27|| for an earlier analysis within 



the color-singlet model). The authors have considered the fragmentation of gluons and 
charm quarks produced at leading order via the Compton- and Bethe-Heitler processes, 
7 + g — > q + g and 7 + g —>■ c+c, respectively. It appears that at large transverse momenta, 
Pt ^ 10 GeV, color-singlet charm quark fragmentation dominates over the photon-gluon 
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Figure 3: Transverse momentum distribution da/dp± for J/tp photo- 
production at the photon-proton centre of mass energy v = 100 GeV. 
The solid line represents the (leading-order) fusion contribution and the 
dashed-dotted line the charm quark fragmentation contribution. The 
dotted and dashed lines represent the gluon fragmentation contributions 
with and without a color-octet component for the S state. The cut on 
the inelasticity parameter z is 0.1 < z < 0.9. From Ref.[26]. 



fusion process while gluon-fragmentation is suppressed by an order of magnitude over 
the whole range of pt, see fig.[3] P6 |. Color-octet contributions to the gluon fragmentation 
function considerably enhance the large- z region, but have no strong effect in the inelastic 
domain z <^ 0.9. Thus, the information that will be obtained from large pt production of 
J ftp particles at HERA can be used to study the charm quark fragmentation mechanism 
and is complementary to the analyses performed at the Tevatron where the large pt 
region is dominated by gluon fragmentation. However, since the cross section for J/tp 
photoproduction at pt ^ 10 GeV is at the most 0(lpb), a large luminosity is required to 
probe fragmentation contributions to quarkonium production at HERA. 



4 Resolved photon contributions 



The photoproduction of a J/ip particle can also take place via a so called resolved photon 
interaction, where the photon couples through one of its hadronic components. For J/ip 
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production within the color-singlet model the process goes like: 

IP -> 9~, + 9 P -»• cc[ 3 S u 1] + g ^ J/ip + g (7) 

meaning that the gluon coming from the photon fuses with the one coming from the proton 
to give a 3 Si color-singlet cc pair which subsequently hadronizes into a J/0. This process 
has been extensively analyzed in the past (see for instance J[|) and found to contribute 
to the overall cross section only marginally everywhere but in the low- z region. More 
contributions, coming from production and radiative decay of x states, are also expected 
to be present. These terms also probe the quark content of the photon. They are expected 
to be of comparable size with the J/ip production process (0) described above (see [[[]] for 
a survey of these and other resolved production mechanisms). 

Within the factorization approach, however, more channels have to be considered, 
where the Jfip production goes via color-octet cc pairs. Namely, the following processes 
contribute: 

9~t + 9 P -> ccpSo.g] +g 
9j + 9 P -^ cc[ 3 S u 8] + g 
9-i + 9 P -^ cc[ 3 Pj,8] +g. 

They are in every respect analogous to the ones which have been argued to greatly increase 
the J/ip production in pp collisions at the Tevatron. They must therefore be carefully 
considered here to understand whether or not they change the picture established by the 
CSM. 

Using the matrix elements squared evaluated in ref.|| and the choice of nonpertur- 
bative parameters displayed in Table |] we can calculate the cross sections and compare 
them with the CSM ones. The results are shown in Table 0, with a minimum p? cut 
needed to screen the collinear singularity present in the 1 Sq and 3 Pj channels. 

These results make clear that the color-octet channels could provide a non-negligible 
increment of the overall J/ip photoproduction cross section. It's therefore worth studying 
in more detail which regions of phase space will be mostly affected. A close look to 
the differential distributions of experimental interest shows that the behaviour of the 
color-octet contributions to the resolved channels is similar to that of the color-singlet 
one. Namely, it is only visible in the low-;? region. Fig.[| shows the z distribution for 
the direct photon color-singlet channel (full line) compared with the old CSM resolved 
photon contribution (dashed line) and the new resolved photon color-octet one (dotted 
line). The resolved contributions can be seen to be enhanced by the color-octet terms, but 
the qualitative picture of them being visible only in the low- z region remains unchanged. 
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a 7 p (nb) 


Channel 




y/s = 100 GeV 








p T > 1 GeV 


Direct 


3 S\. 


1 


13.68 




3 Si. 


1 


.48 




X Sq. 
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.79 






8 


.34 
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1.61 




3p 


8 


.50 




3 p 


8 


2.06 



Table 2: Results for the total cross sections (in nb). 
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Figure 4: Comparison of the inelasticity distributions for direct (singlet) 
and resolved (singlet and octet) J/0 photoproduction. 

5 x photoproduction 

For the production of S"-wave quarkonia, like J/ip, ip' or T, the factorization approach 
coincides with the color-singlet model in the nonrelativistic limit v — > 0. In the case of \ 
bound states, however, color-singlet and color-octet mechanisms contribute at the same 
order in v to annihilation rates and production cross sections and must therefore both be 
included for a consistent calculation pBfl . The production of x states plays a peculiar role 
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in photon-hadron collisions. Unlikely what happens in hadron-hadron collisions, were x s 
are copiously produced by gg, gq and qq interactions, the process 

1 + g^cc[ 3 P J ,l]+9 (8) 

namely the production of a color-singlet cc pair with the x quantum numbers, happens to 
be forbidden at leading order in a s . This is due to the color factor of the two gluons having 
to be symmetric in order to produce a color-singlet states. The gluons could therefore be 
replaced by two colorless photons, and the process \ — > 777 is know to be forbidden by 
the request of charge conjugation invariance. 

To have x production initiated by a direct photon we must therefore either go to 
higher orders within the CSM, adding one gluon to the leading order diagram, or consider 
color-octet mediated channels. 

Within the factorization approach, indeed, x production can still take place at the 
leading 0(aa 2 s ), provided it is a color-octet 3 S\ charm pair which is produced in the hard 
interaction and subsequently hadronizes into a physical x particle. The following two 
processes contribute: 

1 + g^ccfS^+g 

l + q{q)^cc[ 3 S u $ + q{q) 

The first of these two reaction is by far the dominant one in the HERA energy range 
and at low transverse momentum. Within the spirit of the factorization approach, the 
cross section is given by the same short distance cross section evaluated for color-octet 
J/ijj photoproduction (see eq.(|l]) and the second process of eq.(fD), times the appropriate 
matrix elements (O xj [ 3 Si, 8]). These matrix elements have also been fitted to the Tevatron 
data in [7], |[, and found to be of order 10~ 2 GeV 3 . Using, for the sake of simplicity, 
(O X0 [ 3 S u 8]) = KT 2 GeV 3 and the relation 

(0 XJ [ 3 5!,8]) = (2J + l){O X0 [ 3 S 1 ,8}) (9) 

we find 

J2°(lP -> Xj) = 9(O X0 [ 3 S 1 ,8}} x 0.9 nb GeV' 3 ~ 80 pb (10) 
j 

at ^/i =100 GeV and with a minimum px cut of 3 GeV. 

It is worth noticing here a pretty large discrepancy with the result of ref.j29|. Despite 



using a nonperturbative matrix element about a factor of two smaller than ours, it finds 
a {lP ~ * Xi) = 0-13 nb, which is about a factor of five larger than our result. 
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6 Associated J/ip + 7 production 



A particularly distinctive experimental probe of the relevance of color-octet contributions 
in quarkonia production can be the observation of the exclusive process given by the 
associated production of a J/ip and a photon |30|p|: 

ip — > J/i> + 7 (11) 

Within the CSM this process can only undergo via the resolved photon channel, since at 
least two gluons must couple to the heavy quark line to produce a color-singlet cc pair: 

IP -> 2 7 + 9 P -> cc[ 3 S u 1] + 7 -> J/0 + 7 (12) 

A big suppression of the cross section, and the characteristic z distribution typical of 
resolved photon processes, peaked at low z, can therefore be expected. 

Within the factorization approach, on the other hand, other cc states can contribute 
and, in particular, the reaction can now proceed also via a direct photon channel. Indeed, 
the following processes are now possible: 



Direct photon (figj^a): 



1 + 9 P -> cc^S Q ,8] +7 
j + g p ^cc[ 3 S 1 ,8}+ 1 
l + g p ^ cc[ 3 Pj,8] +7 



Resolved photon (fig.[5|b): 



9i + 9p^ cc[ S , 8} +7 

9-y + 9p — > cc[ 3 Si, 1] + 7 (Standard CSM process) 
9j + 9 P ^ cc[ 3 5*1, 8] +7 
^ 7 + ^ P -> cc[ 3 Pj,8] +7 



Though the nonperturbative matrix elements which mediate the hadronization of the 
color-octet cc pairs to a J/ip are suppressed with respect to the color-singlet one (see Table 
[l]), we can however still expect the two following features: 

i) the production of color-octet states via a direct process - rather than a resolved 
one - will at least partially compensate for the smaller matrix elements. We should 
therefore expect an increase in the overall cross section; 

2 see also ref.|pl[ for an analysis within the color-evaporation model. 
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Figure 5: Diagrams contributing within the factorization approach to 
J/ip + photon associated production in direct (a) and resolved (b) photon 
collisions. In (b) the right diagram contributes to 1 5o and 3 Pj states 
production only. 

ii) the z distribution of the J/0, will be more peaked towards one. This is again due 
to the presence of a direct photon coupling as opposed to the resolved one, where 
the g 1 only carries part of the photon energy into the reaction. 

These (and especially the second one) are the reasons why we expect this process to be a 
sensitive probe to color-octet components. 







a 1P (pb) 


Channel 


y/s = 100 GeV 






Vt > 1 GeV 




1 -S'o, 8 




Direct 


3 Si, 8 


7.67 




3 Pj,8 








.35 




3 S 1 ,1 


16.70 




3 S U 8 


.27 


Resolved 


3 Po,S 


1.03 




3 Pi,8 


.14 




3 P 2 ,8 


.97 



Table 3: Results for the total cross sections of J/ip + j photoproduction 
(in pb). At leading-order, the 1 Sq and 3 Pj direct contributions vanish 
identically. 
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The total cross sections (with a minimum cut of 1 GeV) are shown in Table |] for 
7P collisions at a cm energy of 100 GeV (see ref. |3(J for more details). The bulk of the 



cross section can be seen to come from the standard CSM channel gg — > cc[ 3 Si, 1] + 7. 
However, the direct channel gives a non negligible contribution, amounting to about 25% 
of the overall cross section of ~ 27 pb. This increase is however far too small to be 
reliably used to assess the presence of the color-octet terms, given the smallness of this 
cross section (in the picobarn region) and the large theoretical uncertainties involved (like 
the charm mass, the strong coupling and the nonperturbative matrix elements values). 

On the other hand, the study of the differential distributions can make easier to 
disentangle the color-octet contributions from the standard color-singlet one. 

We therefore show in fig.|| the differential distributions related to the total ■yp cross 
sections at y/s = 100 GeV with a minimum p T cut of 1 GeV, presented in Table |^. 
The distributions due to [ 3 Si,I] production in resolved photon collision (continuous line) 
and to [ 3 Si,8] production in direct photon collision (dashed line) only are shown. The 
distributions due to the other color-octet processes do indeed present the same features as 
the ones of [ 3 Si,l], being also produced in resolved photon interactions, but are suppressed 
in magnitude, as can be seen from Table |3]. Their inclusion would therefore not change 
the picture we are going to discuss. 

Fig.| compares the result of the CSM with that due to the production of a color-octet 
3 S\ state in direct photon collision, fig.||a, as predicted by the factorization approach. As 
expected, the effect of the direct photon coupling can easily be seen in at least some of 
the plots. While the pt of the J/ip and the invariant mass distribution of the J/0-7 pair 
are pretty similar for the two processes, the z, rapidity and photon energy distributions 
do indeed show a strikingly different behaviour. 

Recalling that we put ourselves in the so-called "HERA-frame" , with the photon (or 
the electron) traveling in the direction of negative rapidities, we notice how the direct 
photon coupling favours the production of the quarkonium and of the photon in the 
negative rapidities region. This contrasts the case of resolved photon production of color- 
singlet 3 5*i states, which are uniformly produced around the central region. 

As for the z distribution, the resolved photon process predicts a decrease of the cross 
section going towards the high- z region. The direct photon process does on the other hand 
predict the opposite behaviour: the cross section now increases going towards z — 1. The 
small dip in the last few bins is due to the minimum px cut. 

Similarly to the z distribution behaves the photon energy one, which is predicted to 
be much harder in direct photon processes. 

These distributions (which were also checked to be robust with respect to a higher 
cut, so as to be sure of the absence of p™ in effects) could already be good experimental 
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discriminators: observation of a substantial fraction of J/0 + 7 events in the high z region 
would provide good evidence for the presence of color-octet contributions to the overall 
cross section. 
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7 J/ij) production in deep inelastic scattering 



Leptoproduction of quarkonium states has been extensively studied in the framework of 
the color-singlet model 0. Though the color-singlet contribution can explain inelastic 
leptoproduction of J/0 it can not explain the total cross section. Thus, in order to arrive 
at a complete description of J/ip leptoproduction, the color-octet contributions to the 
J/ijj production rate have been calculated in ref. |J3. The authors obtain the following 
expression for the differential subprocess cross section: 



2% 2 e 2 c a s a 2 f dy ( 1 + (1 - yf 



m c s 



,( O «r S0 ,8] ) + < o *i 3p »-a> 3e2+7 ( 2m ^ 



mt 



(O^[ 3 P ,8])8(2m c ) 2 g 2 



mi 



>5(sy-(4m 2 c + Q 2 )) , (13) 



where \/§ is the subprocess center-of-mass energy, Q 2 is the negative invariant mass of 
the photon, and y is the momentum fraction of the J/ip relative to the incoming electron. 
To obtain the cross section one convolutes the expression given in eq. ([13]) with the gluon 
distribution function. The hadronization of the cc pair, produced initially in a color-octet 
state with spin/angular- momentum quantum numbers 2S+1 Lj, into a J/ip bound state 
is parametrized by the NRQCD matrix elements (O Jf *^S , 8]) and (O J ^[ 3 P ,8)}. Note 
that it is precisely the matrix elements appearing in eq. (p~3|) that are at the heart of the 
discrepancy between the CDF measurement and the photoproduction results. 

There is an important point regarding the differential cross section presented in eq. (|T3"D . 
In principle this result is valid for all values of Q 2 , however, there are corrections due to 
higher twist terms that have been neglected in the derivation of the factored form of the 
cross section. These higher twist terms are suppressed by powers of m 2 /Q 2 , and will, 
therefore, vanish for Q 2 3> m 2 . Thus it is necessary to compare the theoretical results 
presented in eq. (p~3"D to experimental data in the large Q 2 regime. In the limit where 
s, Q 2 3> ml, eq.([y|) reduces to 



^ m , Q2 l^ {e9 ^ eJm 



dy 



2-K 2 e 2 c a s a 2 
m c s 



or 

s 



(14) 



Note that for large Q 2 the linear combination (O J ^[ x S , 8]) + 3(O J ^[ 3 P , 8}}/m 2 c is de- 
termined. This is precisely the linear combination of NRQCD matrix elements that is 
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measured at CDF. Therefore J/0 leptoproduction can provide an independent means of 
determining this linear combination of NRQCD matrix elements in a manner different 
from the CDF measurement. 



Currently there exists experimental data on the production of J/0 in fj,N collisions ||33|| , 
however, the values of Q 2 probed in this experiment are too low to be in the asymptotic 
region, and the error on the experimental measurements are too large to allow for an accu- 
rate determination of the color-octet matrix elements. The high-statistics measurements 
to be expected in the future at HERA could definitely help to improve the situation and 
to constrain the color-octet matrix elements. 



8 Conclusion 

We have worked out and reviewed the impact of color-octet contributions on quarkonium 
production in photon-proton collisions and deep-inelastic scattering at HERA. Photo- 
production of J/0, if)', T and x states, associated J/0 + 7 production, fragmentation and 
resolved-photon contributions as well as deep inelastic J/0 production have been discussed. 
We have shown how these reactions can be used to constrain the color-octet matrix ele- 
ments and test the picture of quarkonium production as developed in the context of the 
NRQCD factorization approach. 
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